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; ABSTRACT 

We observed three transits of the extrasolar planet HD 189733b in H I Lyman-a and in a few other lines in the ultraviolet 
with HST/ACS, in the search for atmospheric signatures. We detect a transit signature in the Lyman-a light curve 
with a transit depth of 5.05±0.75%. This depth exceeds the occultation depth produced by the planetary disk alone 
at the 3.5a level (statistical). Other stellar emission lines are less bright, and, taken individually, they do not show the 
transit signature, while the whole spectra redward of the Lyman-Q line has enough photons to show a transit signature 
Q, ' consistent with the absorption by the planetary disk alone. The transit depth's upper limits in the emission lines are 

- ■ 11.1% for Oi A1305A and 5.5% for Cii A1335A lines. 

The presence of an extended exosphere of atomic hydrogen around HD 189733b producing 5% absorption of the full 
. unresolved Lyman-a line flux shows that the planet is losing gas. The Lyman-a light curve is well-fitted by a numerical 

CIh' simulation of escaping hydrogen in which the planetary atoms are pushed by the stellar radiation pressure. We constrain 

' ' the escape rate of atomic hydrogen to be between 10 and 10^^ gs~^ and the ionizing extreme UV fiux between 2 and 

40 times the solar value (l-f), with larger escape rates corresponding to larger EUV flux. The best fit is obtained for 
dM/dt=10^°gs-^ and an EUV fiux fEuv=20 times the solar value. HD 189733b is the second extrasolar planet for 



o 

U , ^ 

^ ■ which atmospheric evaporation has been detected. 
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S^O 1. Introduction lenged by Ben-JafFel (2007), but the apparent discrepancy 

^ . f u I ■ been resolved and the result obtained by Ben- Jaffel 

1.1. Evaporation of hot Jupiters ^^^^ ^^^^^ j^^^ strengthens the evaporation scenario 

' Observations of the extrasolar planet HD209458b in the (Vidal-Madjar et al. 2008). Two subsequent observations 

m Lyman-a line of atomic hydrogen (H i) have revealed that ^t low spectral resolution confirm the evaporation with the 

O , this planet is losing gas (Vidal-Madjar et al. 2003). Indeed, detection of a 5% absorption of the whole Lyman-a hue, as 

O . since the discovery of 51 Peg b (Mayor & Queloz 1995), measured using unresolved emission hue flux during plan- 

the issue of the evaporation of hot Jupiters (also named etary transits observed with the Hubble Space Telescope 

"Pegasides",Guillotetal. 1996) has been raised by Burrows (HST) STIS (Vidal-Madjar et al. 2004) and ACS instru- 

:--] ■ & Lunine (1995) and Guillot et al. (1996). Nonetheless, the ments (Ehrenreich et al. 2008). An independent, new anal- 

/\ ■ discovery of a large number of massive hot JupitersQ and ysis of the low-resolution data set used by Vidal-Madjar et 

b ! very hot JupitertQ has led to the conclusion that the evap- al. (2004) has confirmed that the transit depth in Lyman-a 

oration of massive planets has to be modest. is significantly greater than the transit depth due to the 

In this frame the discovery that the transiting extra- planetary disk alone (Ben- Jaffel, private communication), 

solar planet HD 209458b is indeed losing mass was rather ^11 of these three independent observations show a signifi- 

unexpected. Lyman-a observations of the planetary tran- cant amount of gas at velocities exceeding the planet escape 

sit showed excess absorption due to an extended H i cloud velocity, leading to the conclusion that HD 209458b is evap- 

and provided a lower limit of the H i escape rate on the or- oratmg. 

der of IQio gs-i (Vidal-Madjar et al. 2003; Vidal-Madjar & These observational constraints have been used to de- 

Lecavelier des Etangs 2004). This discovery has been chal- velop several theoretical models to explain and characterize 

the evaporation processes (Lammer et al. 2003; Lecavelier 
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originates from charge exchange between a stehar wind of 
exotic properties and the planetary escaping exosphere (in 
this model the escape is found to be of the same order of 
magnitude), all of these modeling efforts have led to the 
conclusion that most of the EUV and X-ray energy input 
by the host star is used by the atmosphere to escape the 
planetary gravitational potential. Therefore, an "energy di- 
agram" has been developed to describe the evaporation sta- 
tus of extrasolar planets, a diagram in which the potential 
energy of the planet is plotted versus the energy received by 
its upper atmosphere (Lecavelier des Etangs 2007). It has 
also been suggested that, in the case of planets with a mass 
of only a fraction of the mass of a hot Jupiter and orbiting 
close to their parent stars, evaporation can lead to signifi- 
cant modification of the nature of the planet, leading to the 
formation of planetary remnants (Lecavelier des Etangs et 
al. 2004). The recently discovered CoRoT-7b (Leger et al. 
2009) could belong to a new category of planets, which 
was proposed for classification as "chthonian planets" (see 
Lecavelier des Etangs et al. 2004). A review of the evapo- 
ration of hot Jupiters can be found in Ehrenreich (2008). 

However, up to now HD 209458b has remained the only 
planet for which an evaporation process has been observed. 
This raises the question of the evaporation state of other hot 
Jupiters and very hot Jupiters. Is the evaporation specific 
to HD 209458b or general to hot Jupiters? How does the es- 
cape rate depend on the planetary system and stellar char- 
acteristics? The discovery of HD 189733b, a planet tran- 
siting a bright and nearby KO star (V=7.7), offers the un- 
precedented opportunity to answer these questions. Indeed, 
among the stars harboring transiting planets, HD 189733 
belongs to the brightest at Lyman-a; therefore, despite the 
failure of the HST/STIS instrument one year before the 
discovery of HD 189733b, this planet remained one of the 
best targets when searching for atmospheric escape. In this 
paper we present transit observations of HD 189733b made 
with the UV channel of the Advanced Camera for Surveys 
onboard the Hubble Space Telescope (HST/ACS) using the 
same methodology as employed by Ehrenreich et al. (2008) 
for HD 209458b. 

1.2. HD 189733b 

Located just 19.3 parsecs away with a semi-major axis of 
0.03 AU and an orbital period of 2.2 days, HD 189733b 
belongs to the class of "very hot Jupiters" (Bouchy et al. 
2005). More important, since this planet is seen to transit 
its parent star, the planetary transits and eclipses can be 
used to probe the planet's atmosphere and environment 
{e.g., Charbonneau et al. 2008; Desert et al. 2009). 

HD 189733b orbits a bright main-sequence star and 
shows a transit occultation depth of Ri2.4% at opti- 
cal wavelengths (Pont et al. 2007). The planet has 
a mass Mp=1.13 Jupiter masses (Mjup) and a radius 
i?p=1.16 Jupiter radii (i?jup) in the visible (Bakos et al. 
2006; Winn et al. 2007). The short period of the planet 
(2.21858 days) has been measured precisely (Hebrard 
& Lecavelier des Etangs 2006; Knutson et al. 2009). 
Spectropolarimetry has measured the strength and topol- 
ogy of the stellar magnetic field, which reaches up to 
40 G (Moutou et al. 2007). Sodium has been detected 
in the planet's atmosphere by ground-based observations 
(Redfield et al. 2008). Using HST/ACS, Pont et al. (2008) 
detected atmospheric haze, which is interpreted as Mie 



scattering by small particles (Lecavelier des Etangs et al. 
2008b). High signal-to-noise HST/NICMOS observations 
(Sing et al. 2009) have recently shown that the near-IR 
spectrum below 2 /xm is dominated by haze scattering 
rather than by water absorption (Swain et al. 2008). It 
has been tentatively proposed that CO molecules can ex- 
plain the excess absorption seen at 4.5 /im (Charbonneau et 
al. 2008; Desert et al. 2009). Using Spitzer spectroscopy of 
planetary eclipses, the infrared spectra of HD 189733b's at- 
mosphere have revealed signatures of H2 O absorption and 
possibly weather-like variations in the atmospheric condi- 
tions (Grillmair et al. 2008). A sensitive search with GMRT 
has provided very low upper limits on the meter- wavelength 
radio emission from the planet, indicating a weak planetary 
magnetic field (Lecavelier des Etangs et al. 2009). 

2. Observations 

Using the slitless prism spectroscopy mode of the Solar 
Blind Camera (SBC) of the HST/ACS, we observed 3 tran- 
sits of HD 189733b. These observations provide 2D spectral 
images covering the far-UV wavelengths from ~1200A to 
'^1700A. The log of the observations is given in Table[TJ We 
observed during four HST orbits for the first two transits, 
and during three HST orbits for the last transit. For each 
series of exposures within a given HST orbit, we obtained 
four individual exposures with the target at a first position 
on the detector, and four other exposures at a second posi- 
tion (dithering). We thus obtained a total of 32 exposures 
for transits #1 and #2 and 24 exposures for transit #3. 
Each exposure had 257 seconds effective exposure time. 

3. Data analysis 

The data consist of 2D prism dispersed images of the star 
HD 189733. At these far-UV wavelengths, the spectrum of 
a K-type star is composed of well-known bright chromo- 
spheric emission lines arising above a faint photospheric 
continuum. These emission lines form features that are seen 
as bright spots in the 2D images (Fig.[T]). The brightest line 
is the Hi Lyman-a at A1216A. By identifying other features 
(e.g.. Hi, Oi, CII, etc.) the spectrum can be wavelength- 
calibrated. 

The aim of the data analysis is to obtain the best es- 
timate of the brightness of each emission feature (with its 
associated error bar) in order to search for a possible transit 
signature. The first step is to subtract the background of 
each image. This background is caused mainly by the bright 
geocoronal Lyman-a airglow that produces a roughly uni- 
form background entering the detector in the slitless prism 
mode. This background is accurately determined by the 
mean flux measured in each image above and below the 
stellar spectrum. We used two wide zones 100 pixels high 
centered 100 pixels above and below the spectrum. The 
background level is found to vary significantly within a 
given HST orbit. Nonetheless, various tests show that the 
background estimate is robust and its associated error is 
negligible in the final error budget, well below the photon 
noise of all the measured emission features. 

To estimate the brightness of each emission feature, we 
fitted each 2D spectral image with a parametric model and 
calculated the total number of counts in each line after 
subtracting the modeled 2D flux from the other emission 
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Data set 


Date 


Observation 


Transit 






Start 


End 


Start 


End 


Transit #1 


2007-06-10 


02:04:52 


07:35:03 


03:36:59 


05:24:28 


Transit #2 


2007-06-18/19 


22:43:24 


04:12:45 


00:35:58 


02:23:27 


Transit #3 


2008-04-24 


14:01:12 


17:53:22 


15:00:29 


16:47:58 



Table 1. Log of the observations. Time is given in UT. 
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Fig. 2. Plot of the wavelength calibration. The squares 
show the position of each identified emission feature as a 
function of wavelength. The solid line shows the fit by a 
third-degree polynomial. 



Fig. 1. 2D image of the HD 189733 spectrum obtained with 
ACS prism spectroscopy (upper panel). The bottom panel 
shows the 2D-model of this image. In this image, stellar 
emission lines can be seen as bright spots: from left to right, 
the six brightest spots are the Lyman-alpha line, followed 
by the two Oi and Cii emission features, further to the 
right is the Siiv feature, and the two last spots are emission 
by C IV and C i. The stellar continuum emission is seen as 
a horizontal bright line along the spectrum. 
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Fig. 3. Plot of one-dimensional wavelength calibrated ACS 
spectrum of HD 189733. The brightest lines are identified. 



features and from the continuum. For each 2D spectral im- 
age, the fit to the image was performed in 3 steps. First, 
we calculated a preliminary point spread function (PSF) by 
fitting the Lyman-a emission line. In the second step, the 
whole image was fitted with a model of six emission features 
whose positions were free to vary. In this fit, the initial guess 
for the features positions was set close to the six brightest 
spots as seen in Fig. [T] (Hi, Oi, Cii, Siiv, Civ, Ci). This 
second step was used to obtain a preliminary wavelength 
calibration by using the brightest emission lines. Finally, in 



a third step, we performed a fit of the whole 2D image in 
which the list of emission lines was extended to fainter lines. 
In this final stage the model included the wavelengths of the 
emission lines as a fixed parameter, while the wavelength 
calibration was free to vary and was defined by a 3'^'^ order 
polynomial in the X direction and a 1*^* order polynomial in 
the Y direction. The line list input to the model was defined 
by trial and error; few misidentifications at the beginning 
of the process were easily corrected (even for faint lines) 
by checking the lines positions on a wavelength calibration 
plot. We ended up with a line list given in Table [51 
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Wavelength Identification 



(A) 


1 ono n 
izUo.U 


AT T 1 O; TTT 

IN I + bi III 


1 91 7 


H I Ly-Ct 


1240.8 


N V doublet 


1264.7 


Sin* 


1305.4 


Oi* + Oi** 


1335.1 


Cii + Cii* 


1393.7 


Si IV 


1402.8 


Si IV 


1549.5 


C IV doublet 


1656.9 


C I multiplet 



Table 2. List of lines used to fit the 2D images. 



The intrinsic width of the observed emission lines is 
much smaller than the prism resolution, so the lines can be 
used to estimate the instrumental PSF. In all three steps of 
the fitting procedure, the PSF was modeled by the addition 
of six Gaussians. The sharpest Gaussian is an asymmetrical 
2D Gaussian whose widths along the two axes and the ori- 
entation are free to vary. This sharp Gaussian has a width 
of about one pixel, and was found to vary from one expo- 
sure to the next; this is likely due to the telescope jitter. 
The five other Gaussians are axi-symmetrical Gaussians 
of fixed widths : cri =2.5 pixels, cr2=4pixels, (T3=6pixels, 
(T4=10 pixels, and (75=25 pixels. These widths had been ob- 
tained by numerous experiments and allowed us to obtain 
the best fit to the data. 

The stellar continuum was found to be well-fitted by an 
exponential law in the form i^cont(A) = aexp(6A), where 
a and b are free parameters. This stellar continuum model 
was also convolved with the PSF defined in the previous 
paragraph. 

The resulting wavelength calibration (third degree poly- 
nomial) is plotted in Fig. [2] together with the features po- 
sitions obtained at the second step of the fit. It is clear 
that the 3'''^ degree polynomial gives a satisfactory wave- 
length calibration. This calibration also allow correction 
for misidentifications. For instance, we first believed that 
the line between the N v doublet and the O i emission fea- 
ture was a line of Sill from the ground level at A1260.4A. 
But its measured position was significantly off the wave- 
length calibration. We thus realized that this line is in fact 
a line of Sill* at A1264.7A, the Siii line being extincted 
by the interstellar medium absorption. To obtain the wave- 
length calibrated spectrum of HD 189733 plotted in Fig. [31 
we aligned all 2D images using the measured line positions, 
and we calculated the total spectrum by adding the mea- 
sured counts along a band of 5 pixels height. 

With the PSF parameters (8 parameters), the wave- 
length calibration (6 parameters), the brightness of each 
line (10 parameters for a total of 10 lines), the stellar con- 
tinuum (2 parameters), and a free parameter to fit the resid- 
ual background assumed to be constant within the image 
(1 parameter), each image is fitted with a total of 27 pa- 
rameters. 

Finally, after obtaining the best fit to each individual 2D 
image, we calculated a model image for each line using the 
best fit parameters, but changing the brightness parameter 
of the measured line to zero. This model image is subtracted 
from the data image, providing the best estimate of the line 
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Fig. 4. Lyman-a flux as a function of time during three 
transits of HD 189733b. Different colors and symbols are 
for different epochs. The measurements in black diamonds, 
red triangles, and blue squares are for the l'^', 2"^, and 3''^ 
transits, respectively. The curve shows the best fit to the 
data assuming a linear baseline and an optically-thick disk 
transiting the star. The transit depth (defined by the size 
of the occulting disk) is free to vary from one transit to 
another. 



image in which the contribution of all other lines and of the 
continuum has been removed. We measured the brightness 
in the line image using aperture photometry. The error bars 
were calculated using the 2D image of the propagated errors 
provided by the instrument pipeline; these errors are mainly 
dominated by the photon noise. The resulting brightnesses 
were used to obtain the transit light curves plotted and 
discussed in Sect. HI 



4. Transit light curves 

4.1. Hi Lyman-a light curve 

The main objective of these observations is to obtain the 
transit light curve in Hi Lyman-a and to search for an 
atmospheric evaporation signature. For each of the three 
transits, we plotted the resulting light curves (Fig. |4|). We 
fitted each light curve as a function of time by a transit 
model defined by a linear baseline for the stellar Lyman- 
a fiux and the occultation by an optically-thick disk with 
the planetary ephemerides. The results from these three 
independent fits are given in Table [31 

We see that the second transit gives noisier measure- 
ments than the first one, and the observations were made 
only when the planet was off transit or during the partial 
occultation phase. In this transit #2, no measurements dur- 
ing the total phase of the planetary transit were obtained 
because of the Earth occultation of the target as seen by 
HST. In the first and third transits, the data covered all 
phases of the transit light curve, and measurements well 
before and after the transit allowed for estimates of the 
light curve baseline. Measurements obtained during the to- 
tal phase of the planetary transit are well-fitted by a transit 
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Data set 


,,2 
X 


Uegree or 
Freedom 


Absorption 
(%) 


from 0% abs. 


from 2.4% abs. 


Transit #1 


29.4 


29 


7.6 ± 1.5 


5.7 a 


3.8 0- 


Transit #2 


48.8 


29 


5.6 ± 1.7 


3.3 a 


2.0 cr 


Transit #3 


12.7 


21 


2.7 ± 1.2 


2.2 a 


0.2 (7 


All 


97.9 


81 


5.05 ± 0.75 


6.4 a 


3.5 (7 



Table 3. Results from the Lyman-a light curve. 
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Fig. 5. Flux of the Lyman-a line as a function of time dur- 
ing transits of HD 189733b. The flux is normalized by the 
linear baseline as shown in Fig.Sl The curve shows the best 
fit to the data by a single light curve obtained assuming the 
same transit depth for all transits. Symbols are the same 
as in Fig. H) 
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Fig. 6. Same as previous plot in which the data have been 
rebinned by four. The possible difference of transit depth 
between different epochs are visible here. 



light curve of an optically-thick occulting disk (Fig.|31), pro- 
viding an estimate of the size of the H l occulting cloud. 

The full set of Lyman-a measurements can be phase- 
folded to provide a complete transit light curve (Figs. [5] 
and ini). The dispersion of the data rebinned by four is 
roughly divided by two, so there is no sign of red-noise 
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Fig. 7. Plot of the flux in the C ii line for the three HST 
transits. The data are photon-noise limited and the transit 
of the planet is not detected. Only an upper limit on the 
transit depth can be obtained. The solid line shows the 
transit light curve with the 3-cr upper limit on the transit 
depth. 



in the data. The fit of this light curve yields a Lyman-a 
transit depth of 5.05±0.75%. The error bars are estimated 
using the Ax^ method (see e.g., Hebrard et al. 2002), and 
using the same method we can estimate the significance of 
the detection by calculating the difference of between 
the best fit and a fit with no absorption at all (S**^ column 
of TablelSl or a fit with an occultation by the planetary disk 
alone (6*'' column of Table [3]) . Our result shows a transit 
depth that appears to be deeper than the occultation by 
the planetary disk alone (^^2.4%). Using the statistical er- 
ror bars, the excess absorption is at 3.5(7 level. Of course, an 
unknown systematic effect mimicking a transit light curve 
with an excess absorption cannot be excluded; nonetheless, 
the most likely explanation is the presence of a huge cloud 
of atomic hydrogen surrounding HD 189733b. 

The transit depths measured on individual light curves 
of the transits #1 and #3 are different by about 2.5t7. This 
difference is clearly visible in the plot with data rebinned 
by four (Fig. 

4.2. Light curves in other lines 

Using the same procedure as for Lyman-a, we measured the 
brightness of other stellar emission lines and fit the light 
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Line 


A 


N 




Degree of 


Absorption 




3-0" upper limit 




(A) 


measurements 




Freedom 


(%) 


from 2.4% abs. 


(%) 


Ol 


1305 


87 


80.5 


80 


3.7 ± 2.4 


0.5 


11.1 


Cii 


1335 


88 


80.7 


81 


1.0 ± 1.0 


1.3 


5.5 


Siv 


1400 


86 


79.6 


79 


1.6 ± 1.6 


0.8 


11.2 


Civ 


1550 


86 


94.1 


79 


1.7 ± 1.7 


0.5 


7.0 


Ci 


1660 


84 


66.2 


77 


0.7 ± 0.7 


1.5 


6.2 



Table 4. Results from light curves of emission lines other than H i Lyman-a 



curves to measure the transit occultation depths (Table H]). 
In these light curves, few measurements are obviously out- 
liers, in particular in the curve of the C i line. These outliers 
have been removed by flagging the points beyond 3cr of the 
best fit. The number of points used for the fit is given in 
the third column of Table HI 

We detect a transit signature in none of these lines. Only 
upper limit on the transit depths can be derived. For each 
line, the 3-ct upper limit corresponds to the transit depths 
that yield a of 9 above the of the best fit to the data. 
For the two most important species (O i and C ii), we obtain 
a 3-cr upper limit for the transit depth of about 11% and 5%, 
respectively (Fig. [7]). These two species are abundant and 
light ; therefore they should be, after hydrogen, the most 
abundant species in the upper atmosphere. In contrast to 
the case of HD 209458b (Vidal-Madjar et al. 2004), the non- 
detection of oxygen and carbon in the upper atmosphere of 
HD 189733b precludes us from concluding anything about 
the nature of the escape mechanism. Note, however, that 
the lines from the ground level of O i and C ii are expected 
to be strongly absorbed by the interstellar medium (see 
Fig. 2 of Vidal-Madjar et al. 2004), so the observed stel- 
lar emission lines are mainly from the excited levels Oi*, 
Oi**, and Cii*. Therefore, the upper limits quoted above 
should be interpreted with caution, keeping in mind that 
these levels are collisionally populated in gas with typical 
densities above ^lO^cm""^. 



4.3. Total flux redward of Lyman-a 

The planet transit cannot be seen in stellar emission lines 
other than Lyman-a (Sect. H?^ . This is caused by the very 
low number of photons available in those lines. For instance, 
there are about 200 photons per exposure in the O i feature 
at 1305 A; with a dozen of exposures during the planetary 
transit, the accuracy of the transit depth is photon-noise 
limited to about ~2%. This limitation can be overcome by 
the adding photons in a wide zone of the spectra redward 
of the Lyman-a line, including all these other emission lines 
and the stellar continuum together. Using the range from 
1280 A to 1700 A, we find about 3000 photons per exposure 
(corresponding to about 0.5% photon-noise limit on the 
transit depth measurement), to be compared with the 1700 
photons per exposure in the Lyman-a line. We calculated 
the corresponding light curve (Figs. M and ^ and found a 
clear transit signature. In this light curve, we also see a 
flare-like feature superimposed on the transit light curv 



This flare-like feature looks similar to a flare seen at opti- 
cal wavelength during the transit of the planet OGLE- TR- 10b 
(Bentley et al. 2009). This topic will be discussed further in a 
separate paper. 
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Fig. 8. Plot of the flux redward of Lyman-a for the three 
HST transits. Some measurements are clearly outliers; mea- 
surements beyond 3-a from the best fit are plotted with 
stars. Symbols are the same as in Fig. 0] 



This flare took place in the middle of the third visit, at the 
same epochs as outliers already seen in light curves of indi- 
vidual lines, strengthening the need for flagging those mea- 
surements. For the spectra redward of Lyman-a, if outliers 
beyond 3-cr are flagged, we find a transit depth of 2.7±0.6%, 
which is consitent with the transit of the planet with a 
depth of 2.4% (Table [S|). This 4.6cr-detection of the planet 
transit as expected in such data gives confidence to the de- 
tection of the transit seen at Lyman-a. The calculated error 
bars on the depths of both these detected transits are also 
consistent with the theoretical limit from photon noise. 

Finally, the detection of a flare in the middle of the third 
transit could explain why the transit depth in Lyman-a is 
shallower in this last observation. If the measurements ob- 
tained at the epoch of the flare are flagged (the epoch of 
the 3 outliers at 3-(t plus the two adjacent points), we flnd a 
Lyman-a transit depth of 5.7±0.9%. With this last result, 
we would reach exactly the same conclusions. Also, a care- 
ful inspection of individual measurements does not show 
any hint of flare in Lyman-a (Fig. 14]) or in the O i line, in 
contrast to carbon lines (Cii plotted in Fig. [7] and Civ), 
which show suspiciously large flux at the epoch of the flare. 
This behavior is similar to what has been seen for solar 
flares, where large variations are seen in C IV lines, while the 
Lyman-a line that originates in a different part of the stel- 
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Table 5. Results from the light curve redward of Lyman-a 

Redward of Lyman — a 



1.2 r 

- )K 



1.0 




0.9 L , , ^ , , , ^ , , , ^ , , 3 

-2 2 

Time (hours) 

Fig. 9. Plot of the flux redward of Lyman-a. The time has 
been phase-folded and the flux normalized by the linear 
baseline shown in Fig. [S] Outliers at 3-cr are plotted with 
stars. 

lar atmosphere barely shows variations {e.g., Brekket et al. 
1996). We conclude that the flare does not affect the main 
results of the present paper. The detection of the transit of 
the planet in the spectra redward of Lyman-a strengthens 
the results at Lyman-a. 

5. Evaporation of HD 189733b 

5.1. The stellar Lyman-a emission line profile 

As shown in Sect. 14.11 we detect a transit signature in the 
light curve of the H i Lyman-a stellar emission line with a 
depth of ~5%, significantly deeper than the transit 2.4% 
depth due to the planetary disk alone, as seen at optical 
wavelengths. Moreover, the HD 189733b Lyman-a line is a 
broad emission line, with a width of ^±150 km s""'^. From 
the experience with HD 209458b, we know that to inter- 
pret this measurement we must consider the possibility that 
the hydrogen cloud surrounding HD 189733b only absorbs 
a fraction of this stellar line width. Therefore the measured 
5% absorption depth in the total flux of the line is possi- 
bly caused by a larger absorption of only a fraction of the 
line width (see discussion in Vidal-Madjar et al. 2008). To 
interpret the measured absorption depth, it is therefore im- 
portant to have an estimate of the emission line width and 
shape. 

The profile of the Lyman-a stellar emission line seen 
from Earth is composed of the profile of the intrinsic stellar 
line combined with the absorption by the interstellar atomic 
hydrogen and deuterium. Interstellar hydrogen produced a 
narrow absorption at the line center. Deuterium absorption 



2.7 ± 0.6 4.6 a 0.4 a 



50 F 




-200 -100 100 200 

Velocity (km/sec) 

Fig. 10. Plot of the theoretical spectrum of the Lyman- 
a line of HD 189733. The line shape takes the absorption 
by the interstellar medium hydrogen and deuterium into 
account (black thick line). Two different spectra, includ- 
ing absorption by the planetary exophere when it passes 
in front of the star, are superimposed (dashed lines). The 
blue short-dashed line shows the resulting spectrum when 
the escape rate is 3xl0^gs~^ and the ionizing EUV fiux is 
0.1 times the solar value. With these values, in particular 
the low EUV fiux, the blue side of the line is absorbed by 
about 25%, resulting in a decrease of about 13% in the total 
Lyman-a flux. The red long-dashed line shows the theoret- 
ical spectrum of the Lyman-a line during transit for an 
escape rate of 10"'^° gs~^ and an ionizing EUV fiux 20 times 
the solar value. In this last case, the total Lyman-a flux 
decreases by about 5% when the exosphere passes in front 
of the star. 

is weaker and blueshifted by 80kms~"'^. The column density 
of atomic hydrogen toward HD 189733b is not known, but 
a simple relation between the stars' distance and the Hi 
column density can be used. Using the results of Wood et 
al. (2005) and a distance of 19.3 pc for HD 189733b, we infer 
an Hi column of density, A'hi, of about lO^^ '^cm"^. Our 
conclusions are not sensitive to this assumption because 
the column density for a star at ^^20 parsecs is typically 
in the fiat part of the curve of growth. The core of the 
stellar line is fully absorbed by the cold interstellar medium. 
The width of the stellar Lyman-a emission line is obtained 
using the empirical relationship between the star absolute 
magnitude and the chromospheric Lyman-a line FWHM 
(Landsman & Simon 1993). This relation can be written as 
W' = 1.9 X io-^v-/i2.4 .^^hgre the width W- is given in 

A A 

Angstroms. With an absolute V magnitude of 6.1, we infer 
a line width of about 0.61 A, which is FWMH=150 kms~\ 
or for a Gaussian crLyc(=64kms~^. We modeled the double 
peak chromospheric emission line with two Gaussians of 
width ctlj/q/2 and separated by a^ya/i. Combined with 
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the interstellar absorption, we obtained the profile of the 
HD 189733 Lyman-a shown in Fig. \M 



5.2. HD 189733b is evaporating 

We measure a 5% absorption depth over the broad Lyman- 
a line with a width of ~150kms~^, and this can be in- 
terpreted by two extreme scenarios. In the first scenario, 
the hydrogen cloud occults the full line width. In this 
case, the 5% absorption corresponds to a cloud of about 
1.65 Jupiter radii in size with gas velocities reaching over 
150 km s^^. However, at 1.65 Jupiter radii the escape ve- 
locity from HD 189733b is only 49kms~^ (the escape ve- 
locity is 59kms^^ at one planetary radius). In this first 
scenario, the hydrogen gas detected in the Lyman-a tran- 
sit light curve exceeds the escape velocity, and is therefore 
not gravitationally bound to the planet, so the detected gas 
must be escaping. 

In a second scenario, the velocity of the occulting gas 
does not exceed the escape velocity, and so the absorption 
only occurs in the core of the Lyman-a line. Because only 
1/6 of the observed Lyman-a flux comes from the core of 
the line within the 49kms~"'^ radial velocity limit, in this 
scenario the absorption of the core must be at least 30% 
to produce a final absorption of 5% in the full unresolved 
line. A 30% absorption corresponds to a cloud of about 
4 Jupiter radii in size, which significantly exceeds the size 
of the Roche lobe of the planet (3.3 Jupiter radii). In this 
second scenario, the gas detected in the Lyman-a transit 
light curve is beyond the Roche lobe, and is therefore not 
gravitationally bound to the planet, so the detected gas 
must be escaping. 

An intermediate case between these two extreme sce- 
narios can be considered. Even if the Roche lobe was full 
of gas but not beyond its borders, and the gas velocity was 
up to the escape velocity but not above it, the absorption 
depth would only be 3.3%. Therefore, the measured 5% ab- 
sorption depth shows that some gas must be either beyond 
the Roche lobe or at a velocity above the escape velocity. 
The reality is most likely in between the two extreme sce- 
narios with both a fraction of gas outside the Roche lobe 
and a fraction of gas at high velocity. The present obser- 
vations of a 5% absorption of the broad Lyman-a emission 
line demonstrates that HD 189733b must be evaporating. 



5.3. Numerical simulation 

To interpret the Hi light curve in more detail, we mod- 
eled the atmospheric gas escaping HD 189733b with a nu- 
merical simulation including the dynamics of the hydrogen 
atoms. In this N-body simulation, hydrogen atoms are re- 
leased from HD 189733b's upper atmosphere with a random 
initial velocity corresponding to a 10 000 K exobase. This 
temperature corresponds to the temperature expected in 
the upper atmosphere (e.g., Lecavelier des Etangs et al. 
2004; Stone & Proga et al. 2009) and has been actually 
measured for HD 209458b (Ballester et al. 2007). In any 
case, our results do not depend on this assumption because 
atoms are rapidly accelerated by the radiation pressure to 
velocities several times higher than their initial velocities. 
In this model, the gravity of both the star and the planet 
are taken into account. The radiation pressure from the 
stellar Lyman-a emission line is calculated as a function of 




Fig. 11. Plot of the theoretical light curve of the total 
Lyman-a flux for various escape rates and ionizing EUV 
fluxes. The light curve obtained with an escape rate of 
]^QiOgg-i g^j^jj g^j^ ionizing EUV flux of 20 times the solar 
value (thick line) best fits the data (symbols as in Fig. ^ 
with a resulting of 91 for 80 degrees of freedom. The 
dot-dashed line shows the light curve for a planet radius of 
0.1564 times the stellar radius and no additional absorption 
by the planetary atmosphere. 



the radial velocity of the atoms. The self-extinction within 
the cloud is also taken into account. 

The hydrogen atoms are supposed to be ionized by the 
stellar extreme- UV (EUV). In the simulation, their lifetime 
is calculated as a function of the EUV flux, given as an input 
parameter to the model and expressed in solar units. The 
second input parameter is the atomic hydrogen escape rate, 
expressed in grams per second. 

The dynamical model provides a steady state distribu- 
tion of positions and velocities of escaping hydrogen atoms 
in the cloud surrounding HD 189733b. From this informa- 
tion we calculated the corresponding absorption over the 
stellar line and the corresponding transit light curve of the 
total Lyman-a flux to be directly compared with the ob- 
servations. 

5.4. Results 

We calculated a grid of theoretical light curves as a function 
of the two input parameters, namely the escape rate and 
the EUV ionizing flux, and compared them to the measured 
Lyman-a flux. The EUV flux from HD 189733 has never 
been measured and is unknown. However, for a K-type star, 
it can be estimated to be in the range 1 to 100 times the 
solar EUV flux (e.g., Schmitt & Liefke 2004). 

Using the theoretical light curves, we fitted the obser- 
vational data with the escape rate and EUV flux as free 
parameters. For an escape rate below 10'' gs"""^, the oc- 
cultation by the HI cloud is negligible and the resulting 
Lyman-a light curve is similar to the planet occultation 
light curve with a transit depth of ~2.4% (Fig. [TT|) . For a 
large escape rate (dM/dt ^ lO^gs"^) and low EUV flux 
(^3 times solar), the planet has a long trailing Hi cloud 
with a cometary shape, and the occultation is deeper and 
lasts longer than the optical occultation by the planetary 
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Fig. 12. Plot of the absorption depth on the theoretical 
light curve of the stellar H i Lyman-a line as a function of 
the escape rate and EUV ionizing flux. 
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Fig. 13. Plot of error bars for the estimated escape rate and 
EUV ionizing flux. The best fit is obtained for an escape 
rate of lO^'^gs"^ and an ionizing EUV flux of 20 times 
solar. To fit the data at 2cr confidence level, the escape is 
estimated to be between 10^ and 3x 10^^ gs~^ and the EUV 
flux between 5 and 60 times the solar value. The escape rate 
needed to fit the data is correlated with the ionizing flux, 
with large escape rate for a larger ionizing flux. 

disk alone. With these last conditions, the Lyman-a transit 
depth exceeds 10% (Fig.[T2). Both the transit depth and du- 
ration constrains the two parameters. Therefore, the fit to 
the Lyman-a light curve can be used to constrain both the 
escape rate and the ionizing fiux (Fig. [T5)) . We find that the 
data require an escape rate between 10^ and 10^"'^ gs~^ and 
an EUV flux between 10 and 40 times the solar ionizing flux 
(1-cr), with larger escape rates corresponding to larger EUV 
flux (Fig.Hl]). The best fit is obtained for dM/di^lO^" gs'^ 
and EUV flux Feuv=2OF0. With these parameter values, 
the flt to the data with the evaporating planet model yields 
a of 91 for 80 degrees of freedom (Fig. Ilip. This fit is to 
be compared to the fit with a model of a single optically- 
thick disk for which we have a of 97.9 for 81 degrees of 
freedom (Table [3]). Although the number of degrees of free- 
dom is decreased by one in the evaporating planet model 
(because there are two free parameters, whereas only one 



parameter is needed to describe the size of the optically- 
thick disk), this fit is significantly better showing that the 
evaporating planet provides a better scenario to explain the 
observational data. 

6. Discussion 

The model used to interpret the H i transit light curve as- 
sumes a given Lyman-a line width of 0.61 A (Sect 15. ip . 
This value is rather uncertain and might be underestimated 
(Landsman & Simon 1993), and this is strengthened by the 
particularly high chromospheric activity of this star (Henry 
& Winn 2008; Boisse et al. 2009). We tested this assump- 
tion by fitting the data as in Sects. 15.31 and 15.41 but using 
an FWMH of 1.0 A for the simulated HD 189733 Lyman-a 
stellar line. We found that, to fit the observed light curve, 
a broader line requires a lower EUV ionizing fiux, because 
hydrogen atoms need to reach higher velocities to occult 
the broader wings of the line before being ionized. We found 
that an EUV ionizing fiux between 2 and 10 times the solar 
value and similar escape rates as in Sect. [S] provide a good 
fit to the data. Taking the uncertainty on the Lyman-a 
line width into account, we ended up with an EUV ionizing 
fiux between 2 and 40 times the solar value (l-c). This test 
shows that the estimate of the escape rate is not sensitive 
to the assumption on the width of the stellar emission line. 

We constrained the escape rate of atomic hydrogen to 
be about 10^^° gs~^, which is similar to the escape rate pre- 
viously estimated in the case of HD 209458b. The similarity 
between the two escape rates can be explained by compar- 
ing the energy input by the stars and the physical properties 
of the planets. In effect, theoretical models of the evapora- 
tion of hot Jupiters agree to conclude that the escape rate 
can be estimated by assuming that the EUV and X-ray en- 
ergy input from the star is used by the gas to escape the 
gravitation of the planet (Lecavelier des Etangs 2007). The 
potential well of HD 189733b is about twice deeper than 
that of HD 209458b (the planet mass is 1.7 larger, and the 
radius is 12% smaller). However, this well is largely com- 
pensated for by the orbital distance of HD 189733b which is 
1.5 times closer to its star than HD 209458b and the stellar 
type of HD 189733 which as a K-type star is expected to 
provide a higher EUV and X-ray input energy. Within the 
uncertainties, these orders of magnitude show that the sim- 
ilarity between the observed evaporations of HD 189733b 
and HD 209458b results from their similar physical proper- 
ties. In both cases, the observed escape rate is modest, and 
even if the escape of ionized hydrogen is not detectable in 
Lyman-a, evaporation is unlikely to change the nature of 
these two planets: 10^°gs~^ corresponds to the escape of 
about 0.2% of the planets' mass in 10 billion years. 

However, the net effect of variations in EUV flux on the 
observable is not obvious. In effect, although an increased 
stellar EUV flux deposits more energy onto the planetary 
upper atmosphere, and thereby increasing the heating and 
consequently the atmospheric escape (e.g., Lecavelier des 
Etangs et al. 2004), it is also more efficient at ionizing 
the escaping Hi atoms, hence reducing the size of the ob- 
servable Hi cloud measured during the transit. This also 
contributes to explaining why we find comparable absorp- 
tion depths on HD 189733b and HD 209458b, although their 
host stars are of different stellar types. Here it must be re- 
called that the transit depth is measured on the Lyman-a 
line of atomic hydrogen and is interpreted in terms of es- 
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cape rate through a model in which the hydrogen ions are 
produced by ionization by EUV steUar photons; therefore, 
the quoted escape rates do not include escape of hydrogen 
ionized by other mechanisms lower in the atmosphere. 

We did not detect a transit signature of species other 
than hydrogen, because of the photon noise limit on indi- 
vidual lines. The transit signature can be seen by adding all 
photons redward of the Lyman-a including all lines and the 
stellar continuum; in this case the transit depth is consis- 
tent with the absorption by the planetary disk alone. With 
this limitation, it is not possible to conclude anything about 
the escape mechanism, contrary to the case of HD2908458 b 
where detection of oxygen and carbon revealed the hydro- 
dynamical escape ("blow-off") of the atmosphere (Vidal- 
Madjar et al. 2004). 

7. Conclusion 

In summary, we have presented the 3.5cr detection of an ex- 
cess absorption in the Lyman-a transit light curve. This is 
interpreted in terms of the evaporation in the upper atmo- 
sphere of HD 189733b. Any detection at this level of confi- 
dence would call for an observational confirmation, which 
will be possible with the repaired STIS spectrograph. If con- 
firmed, HD 189733b would be the second extrasolar planet 
for which evaporation has been detected. 

In the present observations, the Lyman-a line is not 
resolved, removing any direct information on the radial ve- 
locity of the absorbing gas. This limitation is due to the 
limited spectroscopic capabilities of the ACS instrument 
in the far-ultraviolet, and the ACS was used because of 
the STIS failure. Fortunately, a significant improvement 
is expected after the refurbishment of the Hubble Space 
Telescope performed last May 2009. With the new Cosmic 
Origins Spectrograph (COS) and the repaired STIS, new 
observations should allow us to address the issue of the 
escape mechanism and to probe the dynamics of escaping 
atoms through an improved sensitivity and sampling of the 
absorption along the Lyman-a line profile. 

Finally, the results from each of the three observed tran- 
sits of HD 189733b are not exactly the same (Table[3]). Since 
the observations were obtained at different epochs (the last 
transit was observed ten months later than the first two 
ones), this raises the question of possible time variations 
in the transit depth. From the present data, the detection 
of variations is barely significant ('^2.5cr), and it may not 
be possible to draw definite conclusions on this question. 
Nonetheless, variations would not be surprising, at least 
because of the expected variations in the EUV ionizing 
flux which, for a given escape rate, significantly affects the 
observed transit depth (Fig. [T2]) . Interpretation of obser- 
vations made at a single epoch is always tricky, and any 
estimate or non-detection at a given time should never be 
considered definitive. 
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